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SUMMARY 

Two systems of formulas are presented for the determination of the 
long period perturbations caused by the Sun and the Moon in the motion of 
an artificial satellite. The first system can be used to determine the lunar 
effect for all satellites. The second method is more convenient for finding 
the lunar effect for close satellites and the solar effect for all satellites. 

Knowledge of these effects is essential for determining the stability of 
the satellite orbit. The basic equations of both systems are arranged in a 
form which permits the use of numerical integration. The two theories are 
more accurate and more adaptable to the use of electronic machines than 
the analytical developments obtained previously. 
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LIST OF SYMBOLS 


the mean anomaly at the epoch of the satellite 

the mean anomaly of the satellite for the moment t 

the argument of the perigee of the satellite 

the longitude of the ascending node of the satellite 

the inclination of the satellite orbit with respect to the equator 

the eccentricity of the orbit of the satellite 

the semi-major axis of the orbit of the satellite 

the mean motion of the satellite 

the eccentric anomaly of the satellite 

a (1 - e 2 ) 

the mean anomaly at the epoch of the disturbing body 

the mean anomaly of the disturbing body at the moment t 

the argument of the perigee of the orbit of the disturbing body, 
but with respect to the ecliptic of the fixed epoch 

the longitude of the ascending node of the orbit of the disturbing 
body with respect to ecliptic 

the inclination of the orbit of the disturbing body toward ecliptic 

the eccentricity of the orbit of the disturbing body 

the semi-major axis of the orbit of the disturbing body 

the mean motion of the disturbing body 

the disturbing function 

the position vector of the satellite 

the unit vector directed from the earth’s center to the disturbing 
body 



i, j, k the basic system of unit vectors in the equatorial system of 
coordinate s 


P the unit vector directed from the earth's center toward the 
perigee of the satellite 

R the unit vector standing normally to the orbit plane of the 
satellite in the direction of the angular momentum 

Q the unit vector standing normally to P and R, Q = R x p 

P' the unit vector directed from the earth's center toward the 

perigee of the disturbing body 

R' the unit vector standing normally to the orbit plane of the dis- 
turbing body in the direction of the angular momentum 

Q' a unit vector standing normally to P' and R' Q' = R' x p ' 
h the unit vector along the line of nodes 
F the disturbing force 

F Q the disturbing force F averaged over the revolution of one 
body 

Mq the momentum of F averaged over the revolution of one 
body 

Wq q the momentum of the disturbing force F averaged over the 

revolutions of both bodies 
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ON THE LONG PERIOD LUNI-SOLAR 


EFFECT IN THE MOTION OF AN 
ARTIFICIAL SATELLITE 


by 

Peter Musen 

Goddard Space Flight Center 


INTRODUCTION 


Two systems of formulas are pre sented herein for determining lunar 
and solar long period effects of the first order in the motion of an artificial 
satellite. The first method is based on the theory originally developed by 
Gauss (Reference 1) for a numerical treatment of the very long period ef- 
fects in planetary motion, which was found to be applicable to the case of 
artificial satellites. The second method is valid for close satellites and is 
based on the development of the disturbing function in terms of Legendre 
polynomials . 

Knowledge of these long period effects is essential in determining the 
stability of the orbit and the lifetime of the satellite. As an example, Kozai 
(Reference Z) found that the solar and lunar perturbations have shortened 
the lifetime of Explorer VI ( 1959 Delta) by a factor of ten. The long period 
lunar effect can be deduced by averaging the perturbations with respect to 
the mean anomaly of the Moon and with respect to the mean anomaly of the 
satellite. The first averaging process is performed analytically, the second 
process numerically. The long period solar effect is obtained by averaging 
the perturbations with respect to the mean anomaly of the satellite only. 

Both methods result in the numerical integration of the equations for the 
variation of elements. 
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The interval of integration depends upon the proximity of the satellite to 
the earth and upon the secular changes of the node and perigee produced by 
the earth 1 s oblateness. In a normal case the interval of integration will be 
of the order of several days. For more distant satellites it can be of the 
order of one month or more. The choice of such a large interval is impossi- 
ble if Cowell's method of integration in the rectangular coordinates is used. 

F or more distant satellites it was found that the development of the dis- 
turbing function into series of Legendre polynomials converges so slowly 
that it is impossible to include all important long period terms in an analyt- 
ical development. A large orbital inclination might also contribute substan- 
tially to the slowness of the convergence. In such a case, as for example 
a = 10 earth radii, e = 0.8, the analytical development obtained by Musen and 
Bailie (Reference 3) becomes incomplete. These circumstances gave rise to 
the investigation of the possible use of the Gaussian method for determining 
the lunar long period effects in the motion of an artificial satellite. 

Halphen's form (Reference 4, with corrections by Goriachev, Reference 
5 ) of the Gaussian method was found to be the most convenient , partly because 
it is very adaptable to the use of electronic computers. Some necessary 
modifications were made to include the basic perturbations of the motion of 
the Moon. Also, the Goursat transformation (Reference 6) and the Euler 
summability process were employed to speed up the convergence of the 
hypergeometric series. This method is valid for all values of e, i and a/a\ 

The development in terms of Legendre polynomials can be used if the 
lunar perturbations for close satellites or the solar perturbations for all 
satellites are to be determined. The basic equations of these two cases are 
arranged in the form which permits the use of the numerical integration and 
includes all the significant long period terms. Theoretically, these equations 
are completely equivalent to the extensive analytical developments obtained 
by Kozai (Reference Z) and by Musen and Bailie (Reference 3) but they have 
a more compact and symmetrical form and are more adaptable to the use of 
electronic machines. Thus the formulas given here represent, from a prac- 
tical point of view, a substantial improvement over the previous methods. 


LUNAR DISTURBING FORCE IN THE MOTION OF A 
DISTANT ARTIFICIAL SATELLITE 

Let g 0 , a), n, i, e, a, n be the osculating elements of the artificial 
satellite referred to the earth's equator, and gg oj ' , f i ' f e ', a', n' be 



the elements of the Moon referred to the fixed ecliptic. We assume that 
there is no sharp commensurability between the mean motions n and n ; . It 
is sufficient to take into account the secular changes in the o>' , ft' and to 
neglect the periodic effects and the influence of the precession. Let P be 
the unit vector directed from the center of the earth to the perigee of the 
orbit of the satellite, R be the unit vector standing normally to the orbit 
plane and 


0 = R x p . 

Let P\ Q' , R / be the corresponding vectors of the orbit of the Moon. 
Putting 



”+l 


0 

0 

A x (a) = 

0 


COS CL 

— sin a 


0 

+ 

sin a 

+ cos a 


+ 

cos 

a 

— sin 

a 

0 

+ 

sin 

a 

+ cos 

a 

0 


0 


0 


+ 1 


we have for the components of P, Q, R and P' , Q ; , R' in the equatorial 
system: 

[P, Q, R] - A 3 (H) • A x ( i) • A 3 M , (1) 


[P\ O', R] = A x (e) -A 3 (fl') • A x ( i ' ) -A 3 (a/) . 

The angle between the equator and the ecliptic is e. Designating u 
eccentric anomaly of the satellite, we have for the position vector 


r = Pa(cosu-e) + Q a yl - e 2 sin u . 


U) 

as the 

( 3 ) 


Let 


p = — s + e' P' 

r a 


( 4 ) 


- P s ( COS u 


e) 


Q s yj 1 - e 2 sin u + e 1 P ' 


where s is the parallax, 
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Put 


a. = p - P' , /3 = p • Q' , y = p • R' . 

The system of Halphen's formulas for the imputation of auxiliary quan- 
tities can be slightly modified and rewritten in our notation: 


K 2 = p 2 - 2 + e' 2 , 

K 2 = (1 - e' 2 )(l - a 2 ) - /3 2 - (2 - e ' 2 ) y 2 , 

K 3 = y 2 ( 1 - e' 2 ) 
g 2 = j (Kj 2 - 3K 2 ) , 

g 3 = (2K X 3 ~ 9 K ! K 2 + 27 K 3 ) . 

. _ 27 S 3 


The next step in Halphen's method is the computation of 


and of 






F ( 12 * 12 ' 1,1 ^ ) • 

~ S )’ 


A = 


F 

yr ’ 

i 12 ’ 

5 

144 

— F 

^3 

vr 


9 

g2 3 


tt 82 

\f%2 


11 

12 ' 


2 , 1 


jr r (s) , 




(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 

( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


The values of £ will be proper fractions. The series (Equations Hand 
12) possess an algebraic branch point at if = 0 and the convergence in the 
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neighborhood of this singularity is slow. To obtain more convenient for- 
mulas representing the series rather uniformly in the interval 0 < f < 1 
with the accuracy of at least 10“ 6 , the Goursat transformations and then the 
E summability process are applied. Putting 


a 12 1 h ~ 12 ' 2 “ 1 ^ 


m 


^ a, b, a + b + y, zj 

- i 1 


F(2a,a-b+i, a + b + I, - LzJTZ 


i + v i - r 


we deduce that 


0<f> = 


^3 v i +v/r 


F (i’ b l ' - w ) 


(15) 


Putting 


13 


17 


a 12 ’ b “ 12 


z = 1 - 


F b, a + b - -i , z j 


= ( 1 “Z ) 


we have 


F ( 


F I 2 a “ 1 , a - b + y, a + b 


I izJH. 

2 ' i+^nr 


™ = (rTTfHi- T. *• -)■ 


( 16 ) 


where 


w 


i -_JL 

1 + VF 
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The series (Equations 15 and 16) are alternating and the E-process can be 
applied to speed up their convergence. The forrrula of the E-process for 
hypergeometric series takes the form 


c , - x ) - 


IX 


( a , k ) ( b , k) 

(1, k) (c, k) 


i im y y^ri p 

, .,00 L—k (1, N + j) (c, N + j) Z-J 2 p+ 1 \ J 


( q, k ) = q ( q + 1 ) 


Applying Equation 17 to the series (Equations 15 and 16) and putting N = 3, 
m = 19, we deduce the two following expressions in which the coefficients are 
rapidly decreasing as the power of w increases. 


1 + 


(+ 2.3870942 

- 0.0663082 w 
+ 0.0225632 w 2 

- 0.0117691 w 3 
+ 0.0073743 w 4 

- 0.0051060 w 5 
+ 0.0037250 w 6 

- 0.0027325 w 7 
+ 0.0019070 w 8 

- 0.0011936 w 9 
+ 0.0006337 w 10 

- 0.0002710 w 11 
+ 0.0000884 w 12 

- 0.0000205 w 13 

+ 0.0000030 w 14 

- 0.0000002 w 15 


144 _ . 

4T >A O 


i + 41 ) 

(- 3.7991784 
+ 0 . 3693646 w 

- 0. 1556119 W 2 
+ 0.0889726 w 3 

- 0.0586828 w 4 

+ 0.0419870 w 5 

- 0.0313364 w 6 
+ 0.0233758 w 7 

- 0.0165247 w 8 
+ 0.0104483 w 9 

- 0.0055933 w 10 
+ 0.0024083 w 11 

- 0.0007898 w 12 
+ 0.0001837 w 13 

- 0.0000268 w 14 
+ 0.0000018 w 15 
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The next step is to compute 

K 4 = 9K 3 - K 1 K 2 , 


and 


: 5 = Kj (Kj K 2 - 3 K 3 ) - 2 K 2 2 , 


'll 


K 4 (a 1 ) + + 2 &2 ^3 * 


*22 = K 4 (/3 2 - 1 + e' 2 ) + K 5 + | /2 , 


3 g 2 K 3 


1 - e' 


*3 3 = K 4 y 2 + K 5 + | g 2 [a 2 (l - e' 2 ) + /3 2 - ( 1 - e' 2 )] , 
a 2 1 = K 4 a/3 , 


'12 


*2 3 


*31 


l ll 


*22 


a 32 = (K 4 -f g 2 j/3 y, 
a 1 3 = [k 4 " f g 2 (1 “ «' 2 )] 


ya , 


a 2 - 1 - j Kj , 


P 2 - 1 + e' 2 - I Kj , 


/ _ 
33 

2 1 
y 3 

*1 . 

12 = 

a 2 1 = 

a/3 , 

/ 

23 = 

a 32 = 

Py , 

r 

31 = 

a 1 3 = 

y a . 


Then the matrix is formed with the elements 


A ij = a ij A + a ij B 
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Introducing the dyadic (a matrix), 


O = 


+ A n P'P' + A 12 P'Q' + A 13 P'R' 


+ A 21 Q'P' + A 22 0'Q' + A 23 0'R' 


(18) 


|_ + A 31 RP' + A 32 R'Q' + A 33 R'Rj 

we can represent the "disturbing force," averaged with respect to the mean 
anomaly of the Moon, in the form 


Fo 


2km' 


$ * r 


09 ) 


and the averaged momentum 

w 2km' 

M 0 = 7 T r x O • r . (20) 

a ° 

In addition, in the system adopted in this paper, we have to compute the 
vector (Reference 7) 


K 


o ~ 



( 21 ) 


The vectors and Kq must be averaged along the orbit of the satellite: 


(M 0 0> 

^qq) 



M 0 


d g 



K o d g 


1 

2n 


2jt 





( 22 ) 


(23) 


This averaging is done numerically, by giving the values 0°, 30°, 60° • ■ • , 
330° to u, computing values of M Q (r/a) and K 0 (r/a) for each of these 
angles and then forming the arithmetical means. 
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THE FORM OF VARIATION OF CONSTANTS FOR 
THE LONG PERIOD LUNAR EFFECTS 
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Designating the area "constant" by c , we have 


yfp R . 

(24) 

a ( 1 - e 2 ) , 

(25) 


and 


3 


dc 

dt 


r x F . 


(26) 


Let h be the unit vector along the line of nodes and let (i, j , k) be the unit 
vectors of the equatorial system. Then 


h — i cos H + j sin 0 , 

R = h * k sin i + k cos i 


It follows from Equations 27 and 28 that 


dR 


Uil - i . dfl ^ di 

— - h sin i — + h x R — , 


dt 


and from Equation 26 that 


dt 


dt 


h sin i + h x R ^ 
dt dt 


—E— R + __Lr ib/jL 

1 - e 2 dt dt 


\F 


r x F 


(27) 

(28) 


(29) 


(30) 


The long period effects are deduced from the last equation by applying the 
double process of averaging and taking into account the fact that da/dt = 0: 


di 

dt 


VF 


M, 


00 


h x R , 


(32) 



e 


R • M 00 


de 

dt 


(33) 
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It is of interest to notice that for long period effects, the averaged 
momentum is the only function of the disturbing force contained in the Equa- 
tions 31 through 33. Equations 31 and 32 represent a particular case of 
equations deduced by Makarova (Reference 8) for computing the special 
perturbations of minor planets. From Reference 7, 


d i 
dT 


+ k 


an 

dt 


R ^ 

R dt 




rR • F - 


VF 


RP 


(34) 


where 


k = (* + ?) F -i 7 rr ■ f - <3s) 

We deduce from Equation 34 the long period effect in co , and by taking Equa- 
tion 31 into consideration obtain 

(36) 

whe re 


dco 

d t 


_ \/p p 


K 


00 


cot 1 

Vp" 


M, 


oo 


p ■ k o» = 5? r p ■ K »^ du - 

J o 

O 

Equation 36 is identical in form to the equation for planetary perturbations, 
deduced by Makarova on the basis of different c onsiderations from Musen's 
theory (Reference 7). 

The system (Equations 31, 32, 33 and 36) is to be integrated numerically; 
and the interval of integration can be of the orcer of several days. 


THE LUNAR PERTURBATIONS OF A CLOSE SATELLITE AND 
THE SOLAR PERTURBATIONS OF ALL SATELLITES 


The computation of the lunar perturbations of a close satellite and of 
the solar perturbations for all satellites can b€ accomplished by developing 
the disturbing function into series of Legendre polynomials. For this pur- 
pose the discussion shall be limited to the second Legendre polynomial, and 
the disturbing function takes the form 
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Q 


R 


km * r 2 / 3 9 

r ' 3 \2 cos 



(37) 


where H is the angle between r and r' . Let u° be the unit vector in the 
direction of r'. We have u° with the accuracy up to e' in the periodic 
terms: 


u° - P' (-e' + cos g' + e * cos 2g # ) + Q' (sin g' + e ' sin 2g' ) , (38) 

and Equation 37 becomes 

km 1 


R = -=Z-r ( 3r • u »u° • r - r 2 ) . 


2 -' 3 


(39) 


The "disturbing force" to be used in connection with the variation of 
elements will be 

F - R - ill (3 „V> . , - HO) 

and the momentum of F takes the form 


3km 7 / a^\ 3 
a 73 [r'J 


r x • r , 


(41) 


Substituting 


r - P r cos f + Qr sin f 


into Equation 41, yields: 


M 


3km 7 a 2 


p) 3 [px u°u° • p d a 


+ cos 2 f ) 


+ P x u°u° • Q —~z sin 2 f + Qx U°u° • P ~ sin 


2f 


2 “| 

+ Q * u°u° • Q ^2 ( 1 - cos 2f)J . 


T aking 


( 42 ) 


2rr 

J_ f £ 

2 - I a 


dg = 1 + ^ e 2 , 


_L [ 2v J d 

277 I a 2 C ° 


s 2f dg = + £ 
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we obtain the averaged value of M 0 : 


Mu 


3km' a 2 / ~'^ 3 


2a 


'3 \r 


A.) ( 1 + 4e 2 ) P x u°u° • P + ( 1 - e 2 ) Q x u°u° • Q 


. (43) 


Let <j> , 4> , and 6 be angles which the basic vectors P, Q, and R form with 
the direction to the Moon; hence 


0 = P • U® , cos l/l - 0 • U® , 


e = r • u w . 


(44) 


Replacing M 00 in Equations 31, 32 and 33 by Equation 43, and considering 


we obtain 

. . an 

sin 1 dt 




h 

. p X 

u° = 

+ 

cos 

8 

sin 

u , 



h 

• 0 X 

u° = 

+ 

cos 

8 

cos 

u , 

(h X 

R) 

• 

(P X 

u°) = 

+ 

c os 

8 

cos 

OJ , 

(h X 

R) 

• 

(0 x 

u°) = 

- 

cos 

8 

s in 

OJ , 



R 

• P X 

c 

o 

II 

+ 

cos 

* 

* 




R 

• 0 X 

s 

o 

II 

- 

c os 


t 






h = 

p 

cos 

CO 

i 

O 

sin 


d t 


3km' a 2 

f a ' \ 3 cos 0 

( 1 + 4e 2 ) 

cos 0 sin co + (1-e 2 ) cos i b cos oj 

2a' 3 ' 

\ r / Vp" 



3km' a 2 

/ a' \ 3 cos 8 * 

(1 + 4e 2 ) 

cos <p cos o. “ ( 1 “ e 2 ) cos \p sin a> 

2a' 3 

vw vr 


- 


(45) 


(46) 


^ ^ 2 Ji_ m a - f cos 0 cos . (47) 


dt 


= -Te/^ 


yTT a 


t 3 V r 
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Setting the mass of the earth equal to 1, we have 

k 


= n a , 


and the system (Equations 45 through 47) becomes 


. an 

1 dt 



~~t\ (l + 4e 2 ) cos 0 sin co + (1-e 2 ) cos0coso)jcosl9: 


d i 
d t 



( 1 + 4e 2 ) cos0 cos ^ — (1-e 2 ) cos yp sinco 



In addition we have, with sufficient accuracy, 

' \3 


a ^ - ( 1 - e ' 2 ) 2 + 3e ' cos g' ; 

r ' ' 


and for dco/dt an equation analogous to Equation 36, 


where 


dco Vp" n v . dCl 

dt~ = ‘ ~ P • K o “ cos 1 dF - 


K .if \ll + iW - J-F.rrl 

K 0 2tt J L\ P/ rp J 


dg 


Substituting 


p = tdL (3 u°u° • r - r) 

i- 1 3 


; (48) 


: (49) 


(50) 


(51) 


(52) 
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and with the solar perturbations for all satellites.. This system contains all 
the significant long period terms. In addition to the long period terms, the 
lunar part also contains the terms having periods of the order of one month. 
The interval of integration for the solar perturbations can be taken to be 
approximately one month. However, for the lunar perturbations it would be 
preferable to take the interval to be of the order of days, to obtain a smoother 
curve and to avoid ambiguity in interpreting the results. 


CONCLUSION 


The methods described here are more accurate and more adaptable to 
the use of electronic computers than the analytical development obtained 
previously. These methods can be used in investigations connected with the 
stability of orbits and for the separation of known and unknown long period 
effects in observations. It is recommended that in practical use the analyti- 
cal development be replaced by the semi- analytical solution given here. 
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